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Multivariate density moments (correlators) of arbitrary order are obtained for the multi- 
plicative self-similar cascade. This result is based on the calculation by Greiner, Eggers and 
Lipa (reference M) where the correlators of the logarithms of the particle densities have been 
obtained. The density correlators, more suitable for comparison with multiparticle data, 
appear to have even simpler form than those obtained in Wl. 
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1 Introduction 

It has been shown recently Jl|, |2|, |3| that one can express the multivariate generating function for 
a random multiplicative cascade process in a simple analytical form. Using such a generating 
function one can calculate multivariate density moments (correlators) [f|] of any order and for an 
arbitrary number of bins. It was also argued in |], 0, Q that the cumulant density correlators 
obtained in this way take a particularly simple form. They can be expressed by derivatives of the 
generating function describing one vertex, i.e. a single-step cascade. This beautiful result, however, 
required the generating function and the correlators to be expressed in terms of logarithm of the 
density rather than in terms of the density itself. Although mathematically elegant and adequate, 
e.g., for studies of turbulence, such representation is difficult to use in analysis of multiparticle 
production, as it makes the comparison with the data somewhat complicated. The moments 
expressed directly in terms of the density have been also discussed in || but only the recursive 
formula describing them was given there. 

Indeed, the experimental determinantion of the moments of particle density is greatly sim- 
plified by the observation that they can be approximated by factorial moments of the measured 
particle distribution ||. As this procedure is not applicable for the moments of the logarithm of 
density, the corresponding measurements must involve much more sophisticated methods. This 
difficulty is particularly acute for low multiplicities where Poisson fluctuations of particle number 
can substantially affect the estimate of the logarithm of density. 

In the present note we point out that the generating function derived in fl], ||, §] can also serve 
for calculation of the standard density moments and correlators and thus can be used directly for 
analysis of multiparticle production. The obtained analytic expressions are rather simple: they 
do not even require differentiation but simply evaluation of the generating function at integer 
values of its parameters. They can be considered as a generalization to arbitrary cascades and 
arbitrary correlators of the results obtained in [[| for a particular case of the a-model. Although 
our formulae represent a simple reformulation of those derived in [[l], ^ |3|, we feel that they deserve 
attention because, as explained above, they are better suited for analysis of multiparticle data, 
particularly at moderate multiplicities. 

In the next section we derive directly the formula for multivariate density moments (correlators) 
in a multiplicative cascade. In Section 3 we discuss the examples of one, two and three-point 
correlators. Our conclusions and comments are listed in the last section. 

2 Multi-bin correlators 

We shall use the notation of Ref. [Q. The evolution of the self-similar cascade goes as follows. One 
starts with one bin of a given width A, which is characterised by some quantity e, which can be 



the energy dissipation density, the density of particles etc., depending on the particular process to 
be decsribed by the cascade model. Its value can be chosen to be e = 1 without losing generality. 

This bin, the mother interval, is split into two binsf], the daughter intervals. The contents 
of the two daughter bins is obtained by multiplying that of the mother interval, e = 1, by two 
numbers qo and q\ drawn from the probability distribution p(qo, q\) which is often called splitting 
function and which we consider to be symmetric, p(qo, qi) = p{qi, qo)- The resulting densities are 

eo = 9oe = <7o 

ei = q x e = q x . (1) 

Now the two bins become themselves the mother intervals and are subsequently split into two, 
giving after the second step the densities 

e oo — Q'ooQ'o e oi — QoiQo 

eio = <?io<?i en = 9n9i, ( 2 ) 

where the pairs of multipliers (qoo,qoi) and (qio,qu) are again drawn independently from the 
probability distribution p, the same which was used at the previous step of the cascade. 

After J steps, one obtains 2 J bins, each one addressed by the binary index fci&2 ■ ■ • kj, where 
every kj equals or 1. The density of every bin is the product of J random multipliers q which 
follow the path leading from the original mother interval to this particular bin: 

J 

^k 1 k 2 ...kj — H Qk 1 k 2 ...k J - (3) 

i=i 

Now, the probability distribution of the particular configuration of the densities €k 1 ...k, is a 
product of the splitting functions p(<7fc 1 ...fc-_ 1 o> <7fci..Jfc-ii) taken at each branching point of the 
cascade, convoluted with the appropriate 5 functions: 



p(eo...o, • • • , ei...i) 



j i 

11 \\dqk 1 ...k j - 1 odqk 1 ...k j -iip(qki...k j -iO, fe..ij-ii! 

3=1 fci,...,fcj_i=0 



1 / J 

II *[efci..Jfe J -Il9*i-fcj 

ki,...,kj=0 \ 3=1 



•(4) 



The moments of the logarithms of the densities tk\...kj can be calculated from the generating 
function Z T defined in 0: 

Zt(ct ...o, ■■■, cn...i) = ( exp Yl °k 1 ...k J In e kl ... kj \ , (5) 

\ \ki,...,kj=0 J I 



2 This is the case of the binary cascade. Generalization of all the results obtained here to a cascade characterized 
by three- or more-fold splitting is straightforward 



where the averaging goes over all configurations of the densities occuring with the probability 
Zt provides, by differentiation, the moments K of the logarithms of the densities e: 

K(to...o, ■ ■ ■ ,ti...i) = 

( II ([ln6 fcl ... fcj )]V-A = 

\k 1 ,...,k J =0 I 



(6) 



^...fc^O 



It was shown also in P, [3| that Z? factorizes: 



Zt(&o.. 



C1...1) — 11 II Z T (a kl ...k j _ 1 o, o"fci...fc i _i] 



j=\ fci,...,fei_i=0 



where 



Z T (a , cri) = / dqodqxpiqo, q x ) exp(cr In g + o"i In gi) 
is the binary generating function corresponding to a cascade consisting of only one step and 



j 1 

Vkx...kj = /I z2 a ki...ki 

l=j+l k j+1 ,...,k l =0 



(7) 



(8) 



(9) 



is the sum of all a which are descendant with respect to a kl ,,, kj . For example, (r = 000 + <7 i, 
where <7 o = o"ooo + ^ooi an d o" 01 = a 010 + o" n, and so on till the last step of the cascade. 

Let us now consider the multivariate moments of the densities e kl _ kj themselves: 

M(s ...o, • • • , S1...1) = ( II (tk l ..Je j y h i- h ' ), (10) 

where Sk 1 ...k J can be any non-negative integer numbers^. The average is, again, taken over all 
configurations at the end of the cascade: 



M(so... 0l . . . , S1...1) 



II de kl ... kj {e kl ...k J Y kl - 

k 1 ,...,k J =0 



p(e ...o,---,ei...i), 



:u) 



Comparing (|II]) with fl5|) one notices that taking the generating function Zt, ©, at fixed integer 
points S0...0, • • • ) s i-i gi ves directly the moments M: 

M(s ...o, ■■■, S1...1) = Z T (s Q ...o, . . . , si...i)- (12) 

The factorization (|7f) can thus also be reformulated in terms of the moments M: 



M(s ... , • • • , S1...1) = II II M(s kl ... k] _ lQ} s kl ,.. k] _ xl ). 

j=l ki,...,kj-i=0 



(13) 



3 The formulae below are valid for any non-negative s kl ...kj- Only for integer s kl ._ kj , however, the correlators 
M can be represented by the factorial correlators and thus easily estimated from data. 



The rule (||) holds here similarly for the variables s^..^ ■■ The hierarchy of these variables sitting 
at each link and of the binary functions M corresponding to every branching point is shown in 
Fig. 1 on the example of a 3-step cascade. The binary moments in the r.h.s. of Eq (f[3|) correspond 
to the binary generating function Zt{(Tq,(Ti) taken again at integer points <jq = sq, o\ = Si and 
they read: 

M( So , Sl ) = fdq d qi (qo) so (qi) Sl p(qo,qi)- (14) 



We thus conclude that all the density correlators can be expressed by products of the binary 
moments M(sq, Si). This is surely a strong constraint which can be tested against the experimental 
data on factorial correlators J4|. The simplest possibility of such a test is discussed in the next 



section. 



S 0~ S 00 + S 01 




J 00~ "000 ' "noi 



001 010 °011 100 101 110 



Figure 1: Hierarchy of the integer variables s and the functions M illustrating 
Eq. p 



3 Examples 

Consider first the t -order moment calculated for a single bin addressed by a binary number 
\ = h...lj: 

M(0,0,...,t,...,0) = ((e A ) t ). (15) 

This is the particular case of Eq. ( |T5| ) , where all the variables s K addressed by k = k\ . . . k j 
different from A are set to zero, whereas s\ = t. One can notice that s^..^. = t for any j < J 
and Sk 1 ...k J _ ■ = for all other k\ . . . kj ^ A. Thus, knowing that M(0, 0) = 1, one concludes that 
the moment M(0, . . . , t, . . . , 0) is the product of J binary moments M(t, 0) or M(0, t) which are 
all equal to each other due to symmetry and sit at the vertices of the trajectory leading to the 
considered bin. The resulting moment reads 

(e A )*) = [M(t,0)] J (16) 



and it does not depend on the addres A of the particular bin. Here, we recover the well-known 
result, namely the scaling law of the density moments M. The graphic representation of the result 
(P~6| ) is shown in Fig. 2. 




o o 



J steps 



Figure 2: Graphic representation of the one-point moment of order t for J = 3. The 
moment reads M(0, 0, t, 0, 0, 0, 0, 0) = ((e io)*) = [M(t, 0)] 3 . 



Similarly, one can calculate the two-point correlators of order (£, u) between bins addressed by 
A and /i, 



M(0, 



t. 



,u, 



0) 



fix ^ 



(17) 



The graphic representation of this correlator is shown in Fig. 3 and the result is 

((exYie.y) = [M(t + u,0)} J - d M(t,u)[M(t,0)M(u,0)] d '\ (18) 

where d is the ultrametric distance || between considered bins, i.e. the number of generations 
leading to their common ancestor vertex. In other words, the ultrametric distance between two 
bins addressed by the binary indices A = l\ . . . lj and fi = mi . . . rrij, equals d if 

h = m i f° r J — J ~ d (19) 

and 

h-d+i ^ mj-d+i. (20) 

Calculating the 3-point correlator is now straightforward. Consider the correlator of the order 
(t, u, v) between three bins addressed by indices A, // and u, respectively. It depends only on the 
trajectories leading to the bins. Let us assume that the ultrametric distance between the first 
two bins, d\fj, is smaller than that of the second pair of the bins, d^. This assumption does not 
restrict the generality of the solution. For such ultrametric distances between the bins, the 3-point 
correlator reads: 



fix) [€„) [e t 



[M(t + u + v, 0)} J - d ^M(t + u, v) [M(v, O)]^- 1 [M(t + u, n ^"' 
x M(t,u)[M(t,0)M(u,0)] dx »- 1 . 



fj.v a Xfi 



(21) 




d steps 







Figure 3: Graphic representation of the two-point correlator of order (t, u) for J 
3 and d = 2. The correlator is M(t,0,u, 0,0, 0,0,0) = ((e oo)*( e oio) u ) = M(< 
u,0)M(f,«)M(t,0)M(w,0). 



Eqs (|TE|), (|TB| ) and (|2~I| ) provide, at least in principle, a powerful test of the cascade model. 
Indeed, (|16|) allows to determine the standard intermittency parameter from a linear fit 



In < (e x Y > 



where 



ft 



*-(t) 

lnM(t,0) 



In 2 



(22) 



(23) 



A is the total width of the interval in which the measurement is taken and 5 is the width of the 
bins into which the interval A has been divided. 

Similarly, from (|18|) one obtains the relation: 



1 
m~2 



In 



< (exY >< M M > 



(J - d){f t+u -ft- f u ) + (ft,u -ft- fu) 



where 



Jt,u 



\nM(t,u) 
hi2 



(24) 



(25) 



We see that the r.h.s. of Eq. (^) depends only on J — d, i.e. on the number of steps of the 
cascade leading to the splitting point in which the trajectories of the bins A and /i separate. It 
is thus independent of the bin width S. This feature is very well confirmed by the data 0. The 
reason for this is the fact that the above number of steps is not sensitive to the structure of the 
cascade splittings hidden within the bins of width 5 (this holds provided that the cascade splittings 
conserve the particle density; in the case it does not, one can calculate the appropriate corrections 



Another important feature is that the correlator (^) does not depend on the physical distance 
D between the bins for d = J. It means that for the binary cascade this correlator stays constant 



for distances D larger than A/2, its value being f t)U — f t — f u (For a cascade with tripple splitting 
it would remain constant for D > 2/3 A). It means that knowing the intermittency index f s for 
arbitrary s and the above base- level value of the correlator allows one to determine f ttU . 

The dependence of the correlators on ultrametric distance d is also worth studying. This is not 
easy matter, however, because the relation of d to the distance D between the bins is not simple 
and, moreover, depends on the dimensionality of the problem^. Nevertheless, such an analysis 



seems not out of reach 10 . 


Finally, from Eq. (pl|) one obtains the relation: 


In 2 
(J 
(ft- 


< {e x )\e»ne u y > 


_< (e x y >< {e^Y >< (e v ) v >_ 

~ d^ v )(jt+u+v ~ ft ~ fu — fv) ~ 

-U,V Jt+U Jv) T \Jt,U It 



fu)- (26) 

This relation can allow for checking the consistency of the above results since it does not contain 
any new parameters. As it was noted before, the coefficients f s and f s j can be determined for 
arbitrary s and t using the relations (E3) and 



4 Conclusions 

In conclusion, we have rewritten the results of Ref. |I|, 0] in a form which seems better suited for 
application to data on multiple production. The resulting general formula for density correlators 
of any order can be directly compared to multiparticle data on factorial correlators || |4j] . Since 
such a procedure was already applied to the a model with encouraging results JTTJ] , we feel that it 



may be rewarding to analyze the data in the more general framework of arbitrary cascade model. 
This should allow for a sensitive test of the relevance of the cascade models to the processes of 
particle production and, if they appear to be relevant, for better understanding of the structure 
of the cascade and of its parameters. 

Acknowledgements 

This investigation was supported in part by the KBN Grant No 2 P03B 086 14 and by subsydium 
FNP 1/99. 



This was first explored in Ref. 



References 

[1] M. Greiner, H.C. Eggers, P. Lipa, Phys. Rev. Lett. 80, 5333 (1988). 

[2] M. Greiner, J. Schmiegel, F. Eickemeyer, P. Lipa, H.C. Eggers, Phys. Rev. E58, 554 (1998). 

[3] H.C. Eggers, M. Greiner, P. Lipa, [hep-ph/ 9811204) , proceedings of 8th Internetional Workshop 
on Correlations and Fluctuations CF'98, June 14-21 1998. 

[4] A. Bialas, R. Peschanski, Nucl. Phys. B308, 857 (1988). 

[5] M. Greiner, P. Lipa, P. Carruthers, Phys. Rev. E51, 1948 (1995); M. Greiner, J. Giesemann, 
P. Lipa, P. Carruthers, Z. Phys. C69, 305 (1996). 

[6] A. Bialas, R. Peschanski, Nucl. Phys. B273, 703 (1986). 

[7] V.V. Aivazyan et al. (NA22 collaboration), Phys. Lett. B258, 487 (1991); I. Derado, 
G. Jancso, N. Schmitz, P. Stopa (EMC), Z. Phys. C47, 23 (1990); I. Derado, J. Figiel, 
G. Jansco, N. Schmitz (EMC), Z. Phys. C54, 357 (1992); M.I. Adamovich et al. (EMU-01 
collaboration), Phys. Rev. Lett. 65, 412 (1990); Phys. Lett. B263, 539 (1991); Z. Phys. C49, 
395 (1991); Nucl. Phys. B388, 3 (1992). 



[8] B. Ziaja, |hep-ph/9903444| , [hep-ph/ 9904204 . 



[9] R. Peschanski, J. Seixas, preprint CERN-TH-5903/90 (1990). 

[10] J. Czyzewski, in preparation. 

[11] E. De Wolf, I. Dremin, W. Kittel, Phys. Rep. 270, 1 (1996); M. Charlet, PHD thesis, Nijmegen 
University, The Netherlands 1994. 



